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^ ABSTRACT 
< 

\ The gravitational lensing, as well as the velocity field and the cosmological light-cone 
warp, changes the observed correlation function of high-redshift objects. We present an 

^ analytical expression of 3D correlation function, simultaneously including those three 

^ ' effects. When two objects are separated over several hundreds Mpc along a line of sight, 

0^ ■ the observed correlation function is dominated by the effect of gravitational lensing 

. rather than the intrinsic clustering. For a canonical lambda-CDM model, the lensing 

, signals in the galaxy-galaxy and galaxy-QSO correlations are beyond noise levels in 

i large-scale redshift surveys like the Sloan Digital Sky Survey. 

Subject headings: cosmology: theory — gravitational lensing — large-scale structure of 
O ' universe — methods: statistical 

1. INTRODUCTION 



X 

^ I The correlation function is one of the fundamental quantities in searching into the physical 

origin of the universe. In large-scale redshift surveys, redshifts and spherical positions on the sky of 
luminous objects are used for estimating the spatial distribution of mass, but the former is distorted 
by inhomogeneity of the universe. 

Two intrinsic distortion effect on the correlation function originate in the velocity field (Kaiser 
1987) and the cosmological warp (Ballinger, Peacock & Heavens 1996; Matsubara &: Suto 1996; 
Matsubara 2000). The former comes from the fact that an observed redshift corresponds to the 
recession velocity which is composed not only of the expansion of the universe but also of peculiar 
velocities. The latter distortion is brought about by the nonlinear mapping of the objects from the 
expanding curved space on a light cone onto a flat redshift space (Alcock & Paczyhski 1979). 
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In addition to those intrinsic distortions, there are secondary distortions which are due to 
perturbations of hght rays. While the redshift is altered by the Sachs- Wolfe effect (Sachs & Wolfe 
1967), the spherical position is recast by the gravitational lensing (Schneider et al. 1992; Mellier 
1999; Bartelmann & Schneider 1999). Although the Sachs- Wolfe effect is not so important, the 
gravitational lensing can affect the observable correlation function in forthcoming redshift surveys 
as this effect is efficient for high-redshift objects (Gunn 1967). 

The effect of the gravitational lensing on the angular functions w{9) have been intensively 
investigated so far (Bartelmann & Schneider 1999; Moessner et al. 1998; Kaiser 1992; Villumsen 
1996). Recently, among others, cross correlations of galaxies at different redshifts (Moessner et al. 
1998) are successfully applied to the commissioning data from the Sloan Digital Sky Survey (Jain 
et al. 2000), to single out the weak lensing effect. 

As for the 3D correlation function in redshift space, a qualitative treatment of the lensing effect 
is given in Suto ct al. (1999) to estimate the upper limit of the effect, using the phcnomcnological 
Dyer-Roeder distance (Dyer &; Roeder 1973). It is still not clear whether or not the lensing is actu- 
ally efficient where the intrinsic correlation function is negligible on scales comparable to lOOMpc. 
The main purpose of this letter is to give a quantitative treatment of this issue, consistently includ- 
ing velocity and cosmological distortions, and consequently to show the weak lensing actually has 
detectable effects on 3D correlation function in redshift space when large-scale redshift surveys like 
Sloan Digital Sky Survey (SDSS) are available. 

2. OBSERVABLE QUANTITIES 

We take the homogeneous, isotropic FRW metric with scalar perturbations in longitudinal 
gauge: 

ds^ = a2(r) {-(1 + 2cl))dT^ + (1 - 20) [dx^ + SK\x){d9^ + sin^ ^V)] } , (1) 

where a is the scale factor, and r is the conformal time, dr = dt/a, and Sk is the comoving angular 
distance of the spatial curvature K = fio+Ao — 1- For example, Sk{x) = {—K)~^/'^ sinh [(— i^)^/^x] 
for open universe, K < 0, and Sk{x) = X foi' universe, K = 0. We adopt the unit c = Hq = 1 
throughout this letter. 

From the first-order Einstein equation of the metric in equation (1), the density contrast d{x,T) 
and the velocity field v''{x, r) on scales much less than the curvature scale satisfy 

where H = a/a, f = dlnD/dlna, and D is the linear growth rate (Peebles 1980). The Laplacian 
A is taken with respect to comoving coordinates. 

Let us consider a light ray emitted from an object at comoving coordinates {x,0,(p;t), which 
an observer recieves at tq. The conventional redshift z = a~^ — 1 is given by % = Jq dzH~^, but 
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the actually observed redshift is changed by the line-of-sight peculiar velocity, V = riiv'^, where rii 
is a line-of-sight unit vector, and also by the gravitational potential, (j). From the time-component 
of the geodesic equation of the light ray, the observed redshift Zs is given by (Sachs &: Wolfe 1967) 



= z + {l + z) 



V{x) - V{0) - <P{x) + m - 2 r 

Jro-x 



(3) 



where \vc abbreviate the function on a light cone as V{x) = V{x: d, '^0 ^ x)) and so as <^(x). The 
integral is performed on the light cone for a fixed direction of line of sight. 

Now, wc consider the small angle approximation so that light rays arc confined to a narrow cone 
around the polar axis, 6 <^ 1, and introduce new coordinates 9i = 9cos(p, 02 = 9s'mLp, following 
Kaiser (1998). Adopting the Born approximation, the angular components of the geodesic equation 
reduces to equation for the observed angular components Oga {ci = 1,2): 

esa = 9a + £ dx'SKix - X')da4>{x'), (4) 

where da = SK~^d/d9a- 

The apparent luminosity of the hght is magnified by a factor A{z) = \ dei{d9sh/ d6a)\ = 1 + 2k, 
where k is a local convergence field of weak lensing (Schneider et al. 1992; Bernardeau et al. 1997; 
Kaiser 1998) for a fixed redshift of source object: 

K{z,e) = £dx'9{x,x')dadacl^{x'); g{x,x')^^^^^^^^^^^- (5) 

Due to the magnification, the observed apparent magnitude nig is given by rus = m — 2.51og]^Q A = 
m — 5k/ In 10 (Broadhurst et al. 1995), where m is the apparent magnitude in the absence of lensing. 

The magnitude-limited number density in real space rij-lz, 6; <m) and that in observed redshift 
space ns{zs,9s; <ms) are related by the number conservation equation, n^z'^dzd'^O = UgZg'^dzsd'^Os, 
while the number density is given by n(z, 0; <m) = [1 -|- 5{z, 6)\N{z; <m) / {A-k z^) , where 5{z, 0) is 
the density contrast and N{z; <m) is the magnitude-limited number count per redshift. Evaluating 
the Jacobian, one obtains the relation between n and Us, as well as the relation between density 
contrasts. The result contains the terms with 5, V, (j) and their derivatives. For fluctuations on 
a scale k in units of Hubble distance, such variables scales as V ^ k^^5, dV ^ 5, (f) ^ 
dcf) ~ k~^6, and d'^cp ~ 5, where d schematically represents the spatial derivatives in comoving 
coordinates. Consistently to the small angle approximation, we neglect the fluctuations which 
scale as k"^ and k~^, because k is large enough on scales below Hubble distance. Eventually, the 
distorted density contrast is given by 

S, = 6,-^— + {5a-2)K (6) 

where Sr{z,6) is the number density contrast of the objects in real space, and H{z) is the Hubble 
parameter at z. The logarithmic slope of the number counts a at the limiting magnitude m is given 
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by a{z,m) = d logi^ N {z; <m) / dm (c.f., Moessner, Jain &; Villumsen 1998). The first two terms 
of equation (6) depends strongly on radius z and angles 9, while the last, the lens surface density, 
depends strongly on angles but very weakly on radius. The first term of equation (6) is the real 
density fluctuations, the second term is the velocity distortion (Kaiser 1987; Matsubara & Suto 
1996), and the last term is the effect of weak lensing, which consists of the contribution from the 
modulation by magnification bias, 5aK, and of the alternation of surface density by lensing, —2k. 
In the following, we denote each term in equation (6) as (5s = (5r + (Jy + <^i- 



We consider the correlation function between two objects {zi,6i) and {z2, 02) in the small angle 
approximation, 6 = \6i — ^ 1) and assume zi < Z2 without loss of generality. In the absence of 
the lensing term ^i, the correlation function is given by Matsubara & Suto (1996), which generalize 
the work by Hamilton (1992) to high-redshift objects. In the coordinate system {zi,Z2,9), their 
result, with slight modification allowing difference of the kind of two objects, is expressed as follows: 



where z = {zi + Z2)/2, f5i = f{zi)/bi{zi) ~ fl'^-^{zi)/bi{zi), and bi{zi) is the bias parameter of 
object i [i = 1,2) at redshift Zi. Similarly, a bar for any variable means the evaluation at z 
and subscripts 1,2 assume evaluation at objects 1 and 2, respectively. We denote the comoving 
separation x = \/['S'x(x)]^^^ + (X2 — Xi)^; the comoving cosine /j, = (x2 — Xi)/^^ a^nd P^'s are 
Legendre polynomials, and 



where P{k; z) is the power spectrum at redshift z. This formula is valid only for distant observer 
approximation, x <^ xi ^'^d 9^1. We have omitted the finger-of-God effect which is only 
important on scales less than 10/i^^Mpc or lOOOkm/s. In the following, we are interested in the 
scales of 30Mpc or larger where lensing effect appears, so that we can safely ignore the nonlinear 
effect like finger-of-God effect. 

Adopting the small angle approximation (e.g., Bernardeau et al. 1997; Kaiser 1998; Moessner 
et al. 1998), the correlations involving the lensing term 5\ are obtained as 
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CORRELATION FUNCTION 



Cms = ([<5r(l) + 5v(l)] M2) + 5.{2)]) 

= 1 + l(/3i + /?2) + \piP2\ U^; z)Poifi) 

6 5 J 

- ^(/3i+/92) + ^/3i/?2 Ux;z)P2{i^) + -^PiP2Ux;z)P4{fi), (7) 




(8) 



Cri = (<5r(l)5,(2)) = ^nobi{5a2 - %(X2,Xi)(l +2;i)ep [eSK{xiy,zi] , 



(9) 
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X 



rxi 

/ dxiz)gixi,x)9{X2,x)il + zf^p [9SK{xy,z] , (10) 



Jo 



where dx{z) 



dz/H{z), $,p is the projected correlation function defined by 




Mky)P{k;z), 



(11) 



and $,{r;z) = ^oir',z) is the correlation function in real space at z. The terms ((5i(l)5r(2)), 
((5i(l)(5r(2)) are zero for zi < z^ in the small angle approximation. The velocity-lensing term, 
^vi = (<^v(l)^i(2)) is explicitly calculated to be zero, which is because the term 5^ only depends on 
fluctuations along the line of sight which are smoothed out. If we convolve the above expression of 
^i-l and with a selection function along the line of sight, wc obtain the form of angular correlation 
function with the effect of weak lensing (Bartelmann 1995; Villumscn 1996; Dolag & Bartelmann 
1997; Moessner et al. 1998; Moessner & Jain 1998). 

Thus the total correlation function in 3D redshift space is given by ^tot(-2i, z^-, 6) = ^MS+^ri+6i 
for zi < Z2 and ^ ^ 1. The first term dominates on scales much smaller than Hubble distance, 
while the last two terms dominate on scales comparable to the Hubble distance along the line of 
sight. Therefore, even though ^ms is valid only for zi — Z2 <^ z\ (distant observer approximation), 
we can use the form ^tot even when z\ — zi z\. 

In Figure 1, we plot the total correlation function ^tot together with each component, ^ms, 
^ri, and ^u- With the choice of the CDM-like initial power spectrum (Bardeen et al. 1986) with a 
shape parameter F = 0.25 and a linear amplitude cjg = 1, we use the fitting formula for the fully 
non-linear power spectrum of Peacock &; Dodds (1996) for lensing correlations, ^ri and .^n. Linear 
predictions are also plotted in the lower panels. The nonlinearity of the intrinsic correlation ^ms, 
which is only important for the region z^ — z\'^ 0.003, is ignored. We exemplify the low-density 
fiat model with JIq = 0.3, A = 0.7. In the upper panels, the slope of the number counts is assumed 
as a = 1, and the bias factor as 6 = 1, regardless of the redshift. This example corresponds to 
z ~ 0.2 and m 18 of galaxies as seen in Table 1, in which the slope a is calculated from the 
B-band luminosity function of APM galaxies (Loveday et al. 1992), and of quasars (Boyle et al. 
1988). The limiting magnitudes assumed in Table 1 correspond to estimated SDSS redshift data 
of galaxies and quasars. In practice, the slope a can be observationally determined for individual 
catalogue of quasar or galaxy redshift surveys. Each component of correlations roughly scales as 
^MS Of o'8^bib2, Cii cr8^^o^i(5Q:2 — 2), oc as^^o^{5ai — 2)(5a2 — 2) for other parameters and 
models. 

In the lower panels of Figure 1, the galaxy-galaxy (G-G), galaxy-QSO (G-Q), and QSO-QSO 
(Q-Q) correlations are plotted, assuming the SDSS slope a of Table 1. The bias parameter is set 
6 = 1 and 3 for galaxies and quasars, respectively. The line-of-sight separations are large enough 
in lower panels, so that the intrinsic clustering is negligible. 
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4. DISCUSSION 

The absolute value of the intrinsic clustering component ^ms is a decreasing function of the 
line-of-sight separation, Z2 — zi, except the vicinity of zero crossings. On the other hand, the 
lens-lens component is almost independent on the separation and the density-lens component 
^ri is an increasing function. Those behaviors are understood by the fact that the weak lensing is 
efficient between the object and the observer. Thus, intrinsic clustering component dominates for 
small separations, while lens-lens and/or density-lens components dominates for large separations 
along the direction of line of sight. 

The lower panels in Figure 1 show the region where lensing contribution dominates in the case 
of the SDSS magnitude limits to illustrate the typical magnitude of correlations. Are those lensing 
signals detectable? The statistical uncertainty in estimating the correlation function is given by 
(6$,)^ = n/{SQNiN2) (Peebles 1980), where Ni and N2 are numbers of object in the bin used for 
redshifts zi and Z2, respectively, O is the solid angle subtended by the survey area, and 5Q is the 
fraction in the bin used for angle 6. To increase the signal to noise ratio, it is desirable to use 
large bins for 9. To be specific, wc consider a bin [l',10'], and theoretical curves arc integrated 
accordingly, so that 60, ~ lOOvr [arcmin^]. The effective scale of this bin is given by ^eflf = W /^/2 
for^oc^-^ In the SDSS, n ^ TT [str] ^ 1.27rx 10 [arcmin ], and the estimated numbers of galaxies 
and quasars are Nq = 10^ and Nq = 1.7 x 10^, respectively. Assuming we take sufficiently large 
bins of redshifts (this choice is similar to considering angular correlation functions) , the consequent 
estimates of the statistical error are given by 5.0 x 10~^ for G-G, 8.5 x 10^'^ for G-Q, and 2.9 x 10~^ 
for Q-Q correlations, which arc plotted in lower panels. The S/N ratios turn out to be about 10, 
1.3, and 0.22 for G-G, G-Q, and Q-Q correlations, respectively. 

Therefore, the weak lensing in 3D correlation function of galaxies in the SDSS is definitely 
detectable, and the detection of galaxy-QSO cross-correlation is marginal, while the quasar correla- 
tion by lensing is below the noise level in the SDSS. In order to detect the QSO-QSO lensing effect, 
the sample should be at least 5 times larger than the SDSS, or parameters as, ^o, b, a should be 
larger than assumed values. 

In summary, we have obtained a theoretical prediction of correlation function in redshift space, 
taking into account the effect of weak lensing, together with velocity distortions and cosmological 
distortions on a light-cone. Each effect contributes differently to the correlation function, and is 
realistically detectable. Our result provides a fundamental link between theoretical models and the 
observed correlation function in the 3D redshift survey data. Besides the determination of the power 
spectrum itself, various cosmological parameters, especially the bias parameter, can be estimated 
by proper likelihood analyses, including KL transform of the correlation matrix (Vogeley &; Szalay 
1996; Matsubara et al. 2000). One may also be tempted to assume the cosmological parameters 
before analysing data. In which case, the error originated in choosing the wrong cosmological model 
is roughly given by the order of the redshift z times the error of cosmological parameters, since the 
Alcock-Paczyhski effect is roughly proportional to 2; up to z = 1-2. 
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0[arcmin] 0[arcmin] 0[arcmin] 

Fig. 1. — The correlation function along the line of sight for a flat LCDM model. Long-dashed 
lines: intrinsic clustering, ^mS; dashed lines: lens- lens correlation, ^n, dotted lines: density-lens 
correlation, ^ri, solid lines: total correlation function, ^tot- In the upper panels, zi and are fixed 

in each panel, with a = 1, b = 1. In lower panels, zi and Z2 are fixed and parameters which mimic 
the SDSS redshift catalogue (see text) are assumed. The noise levels for the SDSS are also shown. 
From left to right are plotted galaxy-galaxy, galaxy-QSO, and QSO-QSO correlations. NonUncar 
predictions are plotted except the lines which are not enhanced on small angles in lower panels. 
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Table 1. Slope a{z,m) of the number counts with fixed redshifts. 



Galaxies'^(Siim = 18.8) 



Quasars^(Aim = 20.0) 



z... 
a... 



0.05 
0.20 



0.10 
0.42 



0.15 
0.75 



0.20 
1.2 



0.25 
1.8 



0.30 
2.6 



0.5 1.0 1.5 2.0 2.5 3.0 
0.26 0.29 0.30 0.29 0.29 0.28 



^B-band luminosity function of the APM galaxies (Loveday et al. 1992) is assumed 
^B-band luminosity function of the quasar sample (Boyle et al. 1988) is assumed 



